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Summary 

The traditional fully stressed method performs satisfactorily 
for stress-limited structural design. When this method is 
extended to include displacement limitations in addition to 
stress constraints, it is known as the fully utilized design (FUD). 
Typically, the FUD produces an overdesign, which is the 
primary limitation of this otherwise elegant method. We have 
modified FUD in an attempt to alleviate the limitation. This 
new method, called the modified fully utilized design (MFUD) 
method, has been tested successfully on a number of 
problems that were subjected to multiple loads and had both 
stress and displacement constraints. The solutions obtained 
with MFUD compare favorably with the optimum results that 
can be generated by using nonlinear mathematical program- 
ming techniques. The MFUD method appears to have allevi- 
ated the overdesign condition and offers the simplicity of a 
direct, fully stressed type of design method that is distinctly 
different from optimization and optimality criteria formula- 
tions. The MFUD method is being developed for practicing 
engineers who favor traditional design methods rather than 
methods based on advanced calculus and nonlinear mathemati- 
cal programming techniques. The Integrated Force Method 
(IFM) was found to be the appropriate analysis tool in the 
development of the MFUD method. In this paper, the MFUD 
method and its optimality are examined along with a number of 
illustrative examples. 

Introduction 

The fully stressed design FSD method (ref. 1 ) which is based 
on a simple stress-ratio approach, is an elegant design tool that 
is popular across the civil, mechanical, and aerospace engineer- 
ing industries. However, the FSD is useful only for stress- 
limited designs; it cannot properly handle the displacement 
limitations that have become typical design constraints of 
modem structures. When FSD is extended to handle situations 
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with both stress and displacement constraints, it is called the 
fully utilized design (FUD). Two steps that are required to 
obtain the FUD are (1) generate the FSD for stress constraints 
only, and (2) then uniformly prorate it to obtain the FUD. The 
constant proration factor is obtained to satisfy the single most 
infeasible displacement constraint. Although the FUD thus 
obtained is feasible, it can be an overdesign, which i s the primary 
limitation of the otherwise elegant design method. At present, 
a direct design method to efficiently handle both stress and dis- 
placement constraints is not available. Moreover, sustained 
effort to improve FUD has not been reported in the literature. 
Instead of developing a simpler tool, the designers of the 1 960’ s 
were complicating the approach by applying nonlinear mathe- 
matical programming techniques of operations research 
(refs. 2 to 8) and Langrangian-based optimality criteria meth- 
ods (refs. 9 to 11). Some success has been achieved in design 
optimization; however, these techniques can be computationally 
intensive, and convergence difficulties are frequently encoun- 
tered, even for modest problems (refs. 12 and 13). Despite these 
limitations, design optimization is popular in academia and is 
being improved and promoted for industrial applications, espe- 
cially since there is no alternate design tool that effectively 
handles both stress and displacement constraints. These opti- 
mization methods, to a certain extent, have yet to mature and 
become a standardized design tool for utilization by practicing 
engineers. Imagine the distress of these engineers at finding 
that design has been made more complex by the introduction of 
advanced calculus and variational techniques, without a com- 
parable benefit. Although design optimization is analytically 
elegant, a simpler alternative, such as FSD/FUD, need not be 
abondoned, especially for routine and practical engineering 
design. Further research and development needs to be done 
on direct design methods that do not employ mathematical 
programing techniques. 

This paper outlines the development of a simple FSD/FUD 
type design tool that can handle both stress and displacement 
constraints simultaneously. The proposed design tool is called 
the modified fully utilized design (MFUD) method. In its 
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simplicity, MFUD is comparable to the FUD method, yet it 
alleviates the overdesign limitation that has been associated 
with the traditional FUD technique. When tested on a number 
of problems, MFUD produced solutions comparable to those 
generated by design optimization methods. At this early stage, 
MFUD has been developed for two- and three-dimensional 
trusses with linked member areas as design variables. Success- 
ful completion of the MFUD method for stress and 
displacement constraints will eventually open up an avenue for 
the extension of this method to other types of structures and 
constraints. The proponents of optimization methods can also 
benefit from MFUD by using it to initiate optimization itera- 
tions, thereby alleviating some of the computational burden of 
such methods. 

In this paper the theoretical basis of MFUD is developed and 
illustrated for two examples. A summary of MFUD results, 
along with optimization solutions for several examples, is 
included. The Integrated Force Method (IFM, refs. 14 to 16) is 
shown to be an appropriate analysis tool for deriving the MFUD 
formulas (see appendix A). An analytical examination of the 
optimality of FSD and FUD (see appendix B) is followed by a 
discussion. 


Design Optimization Problem 


Standard nonlinear programming terminology is used to 
formulate the design problem for trusses because solutions 
obtained by the MFUD method are compared with optimiza- 
tion results. The areas of truss members that can be linked for 
practical purposes are considered to be design variables. The 
structures are subjected to multiple load conditions, and con- 
straints are imposed simultaneously on both stresses and dis- 
placements. The number of stress and displacement constraints 
are denoted by J s and J d , respectively, with the total number of 
constraints being m = J s + J d . The J s number of stress constraints 
can be specified as 


Si = 


'JO 


- 1 . 0<0 ; = 1 , 2,. ..,/ 5 


( 1 ) 


where cr. is the stress in the yth member and Oj Q is its permissible 
value. 

Likewise, the J d number of displacement constraints can be 
written as 


SJ s +j ~ 


X J 


1 JO 


-1.0<0 j = 1,2, ...,/<* 


( 2 ) 


where Xj is the yth displacement component and X- is its per- 
missible value. The stress and displacement behavior con- 
straints are feasible provided that g k < 0. 

For a truss with n members, the weight can be considered as 
the objective function for design optimization, and it can be 
written as 


n 

W'({A}) = X^P. A i (3) 

1 = 1 

where p., and A- are the length, density, and area of the ith 
member of the truss, respectively. The computer code auto- 
matically modifies equation (3) for linked design variables, but 
that modification is not elaborated here. 


Fully Utilized Design 

The traditional FUD can be obtained in two steps: (1) gener- 
ation of an FSD and (2) uniform proration of the FSD to obtain 
the FUD. 

An FSD for stress constraints only is generated iteratively by 
using a stress-ratio technique that can be written as 


a? m 1 




i = 1,2 ,...,n 


(4) 


where A^is the area of the ith member at the kt h iteration (unit 
member areas can be used to initiate the iterations). The factor 
R s . for the ith design variable is determined as 

„ max(a u ,a 2i ,...,a L/ ) 

**xfi ~ W 

Oio 

where a Li represents stress in member i for load condition L, 
and a io represents the yield strength of member i. The con- 
verged solution of equation (4) is the FSD, designated as {A Y sd . 
The FSD technique produces very fast convergence, usually in 
about 10 iterations, regardless of problem size. 

Prorating the FSD to satisfy the maximum violated displace- 
ment constraint yields the traditional FUD for simultaneous 
stress and displacement constraints: 




( 6 ) 


where [A ^ is the vector of member areas; g ma x is the value 
of the most violated displacement constraint; and X max and X Q 
are, respectively, the most violated and the allowable displace- 
ment values. The uniform proration factor (1 + g max ) in equa- 
tion (6) produces a feasible design. The FUD is likely to be 
overdesigned because all member areas have been increased by 
the same amount, and it has only one active displacement 
constraint. 

The overdesign condition associated with the traditional 
FUD method can be illustrated by considering displacement 
constraints in the design of a five -bar truss (ref. 17) (see 
Numerical Examples, Example 3). The FUD method pro- 
duces an optimum weight of 62.228 lb, whereas the optimality 
criteria method (OC) and the Sequential Unconstrained Mini- 
mization Technique (SUMT) yield 45.016 and 45.029 lb, 
respectively. In this example, the traditional FUD is 38 percent 
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too heavy. For this problem, the proposed MFUD produced a 
weight of 44.817 lb; this is, respectively, 27.98, 0.44, and 0.47 
percent lighter than the weights produced by the FUD, OC, and 
SUMT methods. In the comparison of MFUD and optimization 
results, more than one optimizer is used because the perfor- 
mance of such methods can be problem-dependent, as is shown 
in reference 12. 


Modified Fully Utilized Design 

The MFUD for simultaneous stress and displacement con- 
straints can be obtained iteratively as follows: 

Step ( 1 ): Identify the design variables to initiate the MFUD 
iterations. The first MFUD iteration can begin from the 
FSD {A \f sd (see eqs. (4) and (5)). For subsequent iterations, the 
areas for stress constraints can be obtained from 

(fA 

A , _ V ! _ W ( = 1,2,..., n (7) 

°io 

where (/\U^ is the maximum force in the fth member for all 
load conditions. This strategy ensures that the final MFUD is 
not biased towards the initial FSD, {A } fsd . 

Step (2): Identify the vq number of violated displacement 

constraints {D} = {g vl ,g v 2 g vq 1 f° r the design obtained in 

step ( 1 ). 

Step (3): Update the design independently for each of the vq 
violated displacement constraints contained in set {D}. (See 
the section Identification of Subset of Design Variables for 
a Violated Displacement Constraint, which shows how only 
a few design variables need to be updated to satisfy a violated 
displacement constraint.) Let the number of design variables 
that should be updated to satisfy a violated displacement 
constraint g vr be qt (where qt < n is the total number of design 
variables). The design update rule for the qith design variable, 
then, is 

Aqj = AfX'i 0 + <?vr ) qi = qlq2,..., qt (8) 

where ^ < 1.0 is a weighted parameter (see the section 
Member Weighted Parameter). The design variables in the 
MFUD method are updated independently, in contrast to uni- 
form proration in case of the traditional FUD method. 

Step (4): Repeat step (3) for all vq numbers of the violated 
displacement constraints to obtain vq design subsets ({A} VI , 
for vi = vl, v2, . . . , vq). 

Step (5): Obtain the design update for the structure as the 
union of the vq designs 

{A} mfud ={A} vl U{A} l ' 2 U...U{A} l " ? (9) 


In the union process, the maximum value should be selected in 
case of member duplication. 

Step (6): Repeat steps (1) to (5) until the design converges. 
The converged design will satisfy both stress and displacement 
constraints. A minimum weight condition is not explicitly 
imposed in the MFUD method; however, as will be shown for 
the examples solved, the weight of the design calculated by the 
MFUD method is very close to the optimum weight generated 
from optimization methods. 

The number of design variables qt that are associated with a 
violated displacement constraint and the weighted parameter 
(see eq. (8)) for each design variable can be easily identi- 
fied when the IFM (refs. 14 to 1 6) is used as the analysis tool. 
The derivation of these two parameters is discussed next. 

Identification of a Subset of Design Variables fora Violated 
Displacement Constraint 

The subset of design variables qt directly associated with a 
violated displacement g vr can be identified by examining the 
displacement-force relationship of IFM (see appendix A): 

{X} = [J][G]{F} (10) 

where {X} is the nodal displacement vector; [/] = m rows of 
[S] -T , with [ S ] being the governing IFM matrix; [G] is the 
diagonal flexibility matrix; and { F) is the internal force vector. 

Because of the banded nature of matrices [/] and [G], for a 
single displacement component X vr equation (10) can be 
expanded as 



The displacement component X vr is an explicit function of the 
subset of member areas referred to by indices q \ , <?2, . . . , qt in 
equation (1 1). Thus it is sufficient to update qt design variables 
to satisfy the violated displacement constraint g vr . Repeating 
this process identifies the critical members for each violated 
displacement constraint in set {/)}. 

Member Weighted Parameter 

Instead of a uniform proration of all the design variables in 
the subset qt for a violated constraint g vr , , individually weighted 
parameters are calculated for each design variable. The deter- 
mination of individually weighted parameters is illustrated by 
the example of a 10-bar truss with a violated displacement 
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constraint at its tip — node 3 along the transverse direction 
shown in figure 1 . In this example, the violated tip displacement 
can be satisfied more effectively by increasing the area of the 
members nearest the support, such as members 1 and 5 in fig- 
ure 1 , rather than the areas of those nearest the free end, such as 
members 3 and 1 0. The weighted parameters bias the design in 
favor of support members 1 and 5. These parameters can be 
selected by examining the sensitivities of the violated displace- 
ment constraints. The sensitivity of violated displacement X vr 
with respect to member area A ^ can be obtained via the IFM 
(ref 19) (see appendix A) as 


a *vr 

3A- 


= J 


vr,qi 


( e,F' 


r \ 

l 

It, 

v 

A vr . 

{ A qi J 


( 12 ) 


For the 10-bar truss example, the displacement sensitivities 
for the four members (i.e., 1,3,5, and 10) are 


^- = - 1101.1 = - 1059.2 

A4, 3 dAj 

( members near the support) 

(13) 

= —67.7 = _336.8 

(members near the free end) 


a violated displacement constraint. The weighted parameter for 
the qith member area can be considered proportional to the 
negative gradient of the violated displacement X vr with respect 
to the member area A*f. 

To safeguard against an overdesign, the weighted parameter 
should also be biased against long members with higher densi- 
ties. Such a condition can be imposed from the gradient of the 
objective function W with respect to member area A •: 

= Pqfiqi 

The gradients of the objective function with respect to 
member areas 1, 3, 5, and 10 (weight density for all members 
is 0.1 lb/in. 3 ) are 


i= 36.0 (%£$ = 36.0 

(members near the support ) 

(15) 

Pi? 3 = 36-0 Ao^io = 

(members near the free end) 

The weighted parameters can be considered inversely propor- 
tional to the gradient of the objective function. 

Equations (12) and (14) can be combined and normalized to 
obtain the weighted parameter for the qith member area: 


Equations (13) show that tip displacement is most sensitive to 
member areas 1 and 5 (near the support) and least sensitive to 
member area 3 (close to the tip). Overall, selecting a proration 
factor proportional to the sensitivities is adequate in satisfying 



Figure 1 . — Ten -bar truss (members are circled, nodes are 
not). 


( \ 

F 


t vr _ 


t 


^ vr , qi 


'V 


max 


( \ 

F 

P E(A"f 


' vr,qm 


' qm 


qm — q\,q2, ...,qt (16) 


Although displacement and weight sensitivities are used in the 
derivation of equation (16), their calculations require trivial 
computational effort since the force vector {F} and displace- 
ment coefficient matrix [7] are available from analysis. 

The final normalized weighted parameters for member areas 
1, 3, 5, and 10 of the 10-bar truss are 


£ 3 = 1.000 ^5 = 0.962 

(members near the support) 


= 0.059 £, 0 = 0.216 

(members near the free end) 


( 17 ) 
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The area for the qith member of the vrth violated displace- 
ment constraint can be updated by using the following formula: 

A V = A qfz>q r i t p{' + 8vr) qi = qhq2,...,ql (18) 

The preceding equation is identical to equation (8) except for a 
step length t p < 1 . The step length restricts member areas against 
rapid changes, for which the current analysis may no longer be 
valid. Default step lengths of t = 0.5 for constraint violations 
exceeding unity (i.e., g vr > 1.0) and t p = 0.25g vr for other 
constraint violations were found satisfactory for most problems 
that have been solved. 

The design obtained from equation ( 1 8) for the displacement 
component X vr is designated by the vector {A} vr . Likewise, 
areas can be obtained for all vq of the violated displacement 
constraints contained in set {Z>}. The union of the vq area 
subsets produces the area vector {A r for the truss. This 
area is then compared with the area for stress constraints 
(eq. (7)), and the final MFUD area (A} m -^fora&th iteration is 
obtained as 

Af ud = ma x[A{ sd y A^ x ] I = 1.2.. ... n (19) 

The MFUD iterations are continued until convergence is 
achieved for both stress and displacement constraints. 

Numerical Examples 

The MFUD method was applied to a number of examples. 
The solutions are compared with those obtained by the tradi- 
tional FTJD method and with the optimum solutions generated 
by several optimization algorithms, such as SUMT (ref. 20); 
Sequential Quadratic Programming, SQPJMSL (ref. 21); FD 
(ref. 22); and OC (ref. 9). The MFUD process is illustrated in 
detail here for the first two examples, a three-bar truss and a 
cantilevered truss, under a wide range of linked displacement 
constraints. For other examples, only the final results (summa- 
rized) are presented. 

Example 1: Three-Bar Truss 

A three-bar truss with Young’s modulus E = 30 000 ksi, 
density = 0.10 lb/in. 3 , and allowable strength a o - 20 ksi is 
depicted in figure 2. The truss is subjected to two load condi- 
tions; the first has two load components (P x = -50 kips and P y 
= -100 kips), whereas the second has only one component (P x 
= 50.0 kips). The truss has 10 behavior constraints, consisting 
of 3 stress and 2 displacement constraints (at node 1, 
X ljC <0.2 in. and X ly <0.05 in.) for each loadcase. The optimum 
solution for the three-bar truss was generated by using three 


« — 100 in. — — 100 in. — ►] 
2 3 4 



Figure 2. — Three-bar truss (members are 
circled, nodes are not). 


optimizers: SUMT, FD, and OC. Initial designs of unity were 
used for all design methods. The SUMT and FD optimizers 
converged to optimum weights of 100.07 and 99.95 lb, respec- 
tively, whereas the OC optimizer generated a slight overdesign, 
reflected in a weight of 101.33 lb. The FSD (for stress con- 
straints only), as determined by the stress-ratio technique, gave 
Aj = 1.182, A 2 = 2.504, and A 3 = 3,533 in. 2 The FSD violated 
one displacement constraint (X ly ) under the first load condi- 
tion. The traditional FUD, which satisfied the violated con- 
straint, gave Aj = 1.574, A 2 = 3.336, and A 3 = 4.706 in 2 The 
FUD had only one active displacement constraint and was 
overdesigned by 22.2 percent, with a weight of 122.182 lb. 

The MFUD for the truss converged to an optimum weight of 
99.97 lb (see table I). The MFUD results compare well with 
those generated by SUMT, FD, and OC optimizers (see 
table II). The MFUD, SUMT, and FD methods yielded identi- 
cal numbers of active stress and displacement constraints; 
however, the OC method produced only one active stress 
constraint (a one-fifth of 1 -percent constraint thickness is 
considered active). Overall, the MFUD method performed 
satisfactorily for this problem. 

The convergence characteristics of MFUD, along with those 
for SUMT, FD, and OC are depicted in figure 3. MFUD 
converged rapidly and monotonically in 24 reanalysis cycles 
that included the 12 reanalyses to obtain the FSD. The conver- 
gence characteristics for FD were rather uneven, requiring 47 
reanalysis cycles to reach the optimum solution. SUMT andOC 
solutions required 62 and 80 reanalysis cycles, respectively. 


TABLE I.— THREE-BAR TRUSS RESULTS FOR A FEW 
MFUD ITERATIONS 


Iterations 

Weight, 

lb 

Member area. 

in. 2 

Violated 

constraint. 

A i 


A i 

a 0 

91.714 

1.182 

2.502 

3.533 

0.3322 

1 

93.770 

1.179 

2.710 

3.534 

.2626 

2 

94.669 

1.149 

3.042 

3.445 

2131 

10 

99.000 

1.100 

3.673 

3.330 

.0329 

14 (final) 

99.966 

1.088 

3.841 

3.265 

.0009 


“Represents FSD. 
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TABLE II. — THREE-BAR TRUSS: RELATIONSHIP OF ACTIVE STRESS AND 


DISPLACEMENT CONSTRAINTS TO OPTIMUM DESIGN 


Design 

method 

Weight, 

lb 

Member area, in. 2 

Number of active 
constraints 

Number of 
reanalyses 

A, 

A, 

A 3 

Stress 

Displacement 

MFUD 

99.97 

1.088 

3.841 

3.265 

2 

1 

24 

FUD 

122.18 

1.574 

3.336 

4.706 

- 

1 

10 

SUMT 

100.07 

1.088 

3.848 

3.267 

2 

1 

62 

FD 

99.95 

1.092 

3.855 

3.250 

2 

1 

47 

OC 

101.33 

1.053 

3.913 

3.345 

1 

- 

80 



Figure 3. — Convergence history of three-bar truss. 


Example 2: Cantilevered Truss With Stress and a Sequence 
of Displacement Constraints 

The 10-bay cantilevered steel truss shown in figure 4 is the 
second example. The truss has 22 nodes and 50 members. It is 
made of a material with a Young’s modulus E = 30.0x1 0 6 psi, 
a weight density p = 0.289 lb/in. 3 , and an allowable strength 
a o = 2Gxl0 3 psi. The truss is subjected to two load conditions. 
The first is a single load of 100 kips at node 22 along the nega- 
tive y-direction. The second is a load of 10 kips applied at all 
nodes along the top chord (i.e., node numbers 4, 6, ... , 22; see 
fig. 4). Stress constraints are considered for all 50 bars — 100 
stress constraints for both load conditions. Displacements along 
the transverse direction are specified at the free end (nodes 21 


and 22) as well as at the center span (node 1 l).The magnitude of 
the center span displacement is specified to be a quarter of the 
tip displacement. The tip displacement is a parameter that 
changes, ranging in magnitude between 0.05 and 1 .5 in. The 
stress constraints dominate the design only when the magni- 
tude of the tip displacement exceeds 1.4 in. When displace- 
ments are less than 1 .4 in., both stress and displacement become 
active constraints. Constraint activities of final designs ob- 
tained with MFUD, FUD, FD, and SQPJMSL are given in 
table III for three design situations (tip displacements of 0.5, 
1 .0, and 1 .5 in.). 

Optimum weights obtained by MFUD, FUD, and FD methods 
were normalized with respect to the weight obtained by the 
SQPJMSL method, which performed best for the entire dis- 
placement range 0. 05 <5< 1.5 (see fig. 5). Notice the following 
observations from figure 5 and table III. The MFUD method 
performs adequately for the entire displacement range, with the 
error not exceeding 2 percent of the optimum solution gener- 
ated by SQPJMSL. In contrast, the maximum error in the weight 
obtained by the FUD method was about 26 percent. When 
stresses dominated the design (i.e., displacement values 
exceeded 1 .4 in.), the FD optimizer produced a 4-percent over- 
design, but MFUD, FUD, and SQPJMSL converged to the 
correct optimum. In other words, the traditional FUD method 
exhibited an overdesign condition when the displacements were 
active constraints, but converged to the correct solutions when 
the displacements became passive constraints, that is, at the FSD 
condition. Optimizers SQPJMSL and FD, and the MFUD 
method provided the same number of active displacement con- 
straints for the entire displacement range (see table III). The 
MFUD method, however, produced a greater number of active 



10 members at 20 in./member = 200 in, 


Figure 4. — Cantilever truss (members are circled, nodes are not). Displacement limitation 
S at center is one-fourth of that at the end. 
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TABLE III — CANTILEVERED TRUSS: 
RELATIONSHIP OF ACTIVE CONSTRAINTS TO 
OPTIMUM WEIGHT OF TRUSS AT VARIOUS 
DISPLACEMENTS 


Design method 

Number of active 
constraints 

Weight, 

klb 

Stress | Displacement 


At d 

= 0.5 in. a 


MFUD 

4 

2 

4.701 

FUD 

- 

1 

5.855 

FD 

- 

2 

4.669 

SQPJMSL 

2 

2 

4.618 


At d 

= 1 .0 in. 


MFUD 

23 

1 

2.521 

FUD 

- 

1 

2.928 

FD 

20 

1 

2.540 

SQPJMSL 

21 

1 

2.538 


At d : 

= 1.5 in. 


MFUD 

44 

- 

2.196 

FUD 

44 

- 

2.196 

FD 

25 

- 

2.279 

SQPJMSL 

41 

- 

2.191 


“d = tip displacement. 


SQP-IMSL 



Figure 5. — Weight normalized with respect to SQP-IMSL 
as a function of displacement constraints of a cantilever 
truss 


stress constraints. For example, when the displacement limita- 
tion was 0.5 in., the number of stress constraints produced by 
MFUD, SQPJMSL, and FD were 4, 2, and 0, respectively. For 
this problem, the MFUD method produced results equal or 
superior to the FUD method for all three displacement limita- 
tions (0.5, 1 .0, and 1 .5 in.), as shown in table m. A comparison 
of the weights obtained by the MFUD method and the 
SQPJMSL optimizer shows that SQPJMSL outperforms 
MFUD by 1.8 percent at a displacement of 0.5 in. and by 
0.2 percent at a displacement of 1 .5 in. When the displacement 
value is 1 .0 in., the MFUD weight is better than the SQPJMSL 
weight by 0.69 percent. 


Example 3: Five-Bar Truss 

The five-bar aluminum truss (refs. 17 and 18) depicted in 
figure 6 was subjected to a single load and had a single displace- 
ment constraint in the transverse direction at node 4. The design 
parameters obtained by the four methods are summarized in 
table IV. For this truss, MFUD produced results slightly superior 
to SUMT and OC. The traditional FUD produced a design that 
was 39 percent too heavy. 



Figure 6. — Five-bar truss (members are circled, nodes are 
not). 


TABLE IV.— FIVE-BAR TRUSS DESIGN 
RESULTS 


Results 

MFUD 

FUD 

SUMT 

OC 

Optimum weight, lb 

44.817 

62.228 

45.029 

45.016 

Member area, in. 2 






A, 

0.001 

1.068 

0.001 

0.001 

a 2 

1.475 



1.501 

1.499 

A 3 

0.001 



0.001 

0.001 

a. 

2.124 



2.119 

2.120 

A, 

0.001 


r 

0.001 

0.001 

Active displacement 






constraints 

1 



1 

1 


Example 4: Tapered Five-Bar Truss 

The tapered five-bar steel truss (ref. 23) shown in figure 7 
was subjected to two load conditions and had five stress and two 
displacement constraints for each load condition. The attri butes 
of the designs generated by the four methods are summarized 
in table V; all the optimum weights are in good agreement. The 
MFUD method produced five active constraints, whereas SUMT 
and OC produced three and four active constraints, respectively. 


NASA TM— 4743 


7 



♦ v 



Figure 7. — Tapered five-bar truss (members are circled, 
nodes are not). 


TABLE V.— TAPERED FIVE-BAR TRUSS DESIGN 
RESULTS 


Results 

MFUD 

FUD 

SUMT 

OC 

Optimum weight, lb 

6528.72 

6549.67 

6541.52 

6549.02 

Member area, in. 2 






27.16 

27.38 

26.23 

26.47 

a 2 

9.37 

9.42 

10.37 

10.33 

a 3 

22.04 

21.91 

21.35 

21.23 

a 4 

11.18 

11.27 

11.91 

11.97 

a 5 

1.63 

1.63 

1.81 

1.46 

Active constraints 





Stress 

3 

— 

2 

3 

Displacement 

2 

1 

1 

1 


Example 5: Ten-Bar Truss 

The popular 10-bar truss (ref. 17) depicted in figure 1 was 
subjected to a single load and had 1 0 stress and 4 displacement 
constraints. Table VI summarizes the design results 
obtained by the four methods. MFUD produced a weight 
2 percent higher than that generated by the SUMT optimizer. 
The design generated by the traditional FUD, however, was 
13 percent heavier than the SUMT optimum design weight. 


TABLE VI — TEN-BAR TRUSS DESIGN RESULTS 


Results 

MFUD 

FUD 

SUMT 

OC : 

Optimum weight, lb 
Member area, in. 2 

5164.11 

5741.21 

5057.51 

5061.86 

A, 

22.97 

29.11 

23.19 

23.54 

a 3 

0.33 

0.36 

0.55 

053 

A, 

31.50 

28.58 

30.46 

30.86 

A 8 

7.43 

20.78 

7.43 

7.48 

Aio 

Active constraints 

21.58 

20.03 

21.64 

21.09 

Stress 

— 

— 

1 

— 

Displacement 

2 

1 

2 

2 


Example 6: Tapered Ten -Bar Truss 

The tapered 10-bar aluminum truss (ref. 23) depicted in 
figure 8 was subjected to two load conditions, each with 10 
stress and 4 displacement constraints. Table VII presents the 
design results produced by the four different methods. For this 
example, SQP_IMSL provided the best optimum weight; SUMT 
and MFUD designs were 0. 1 and 0.4 percent heavier, respec- 
tively. The active constraints for MFUD, SUMT, and 
SQP_IMSL numbered 5, 7, and 7, respectively. 





Figure 8. — Tapered ten-truss (members are circled, 
nodes are not). 


TABLE VII.— TAPERED TEN-BAR TRUSS DESIGN RESULTS 


Results 

MFUD 

FUD 

SUMT 

SQPJMSL 

Optimum weight, lb 

3272.64 

3350.60 

3260.75 

3258.26 

Member area, in. 2 





A, 

58.55 

62.49 

55.97 

54.91 

a 3 

2.29 

1.95 

1.98 

2.36 

a 5 

34.87 

35.73 

38.25 

40.07 

a 8 

19.67 

20.11 

21.15 

22.28 

Ajo 

5.84 

7.27 

6.78 

529 

Active constraints 





Stress 

3 

— 

5 

5 

Displacement 

2 

1 

2 

2 


Example 7: Twenty-Five-Bar Truss 

The 25-bar aluminum truss (ref. 7) in figure 9 was subjected 
to 2 load conditions, with 25 stress and 6 displacement con- 
straints for each load case. The bars’ areas were linked to obtain 
8 independent design variables. The attributes for the optimum 
design for this truss are summarized in table VIII. SUMT, FD, 
and MFUD produced comparable optimum weights; however, 
the MFUD weight was 0.23 percent lighter than that of the FD 
optimizer. The active constraints for MFUD, SUMT, and FD 
were 8, 6, and 8, respectively. The weight generated by the 
traditional FUD method was 6.4 percent heavier than that of the 
FD optimizer with a single active constraint. 
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TABLE VIII.— TWENTY-FIVE-BAR TRUSS 
DESIGN RESULTS 


Results 

MFUD 

FUD 

SUMT 

FD 

Optimum weight, lb 
Member area, in . 2 

380.26 

404.44 

381.71 

381.12 

A, 

0.01 

0.02 

0.01 

0.01 

A, 

2.07 

2.27 

2.06 

2.11 

A, 

0.01 

0.01 

0.01 

0.01 

A, 

1.16 

1.34 

1.16 

1.17 

A„ 

Active constraints 

1.86 

1.72 

1.88 

1.83 

Stress 

4 

— 

2 

4 

Displacement 

4 

1 

4 

4 


Figure 9. — Twenty five-bar truss (members are circled, nodes 
are not). 


Example 8: Simply Supported Truss 

Figure 10 shows a 10-bay steel truss with 51 members sub- 
jected to a single load. All bar areas were considered indepen- 
dent variables. The results obtained for 5 1 stress and 2 midspan 
transverse displacement constraints are summarized in table IX. 
For this example, the MFTJD weight lies between the optimum 
weights generated by the FD and SUMT optimizers. 


2 [T0-A-^OE^ 


40 in. 



10 members at 20 inVmember = 200 in. 


Figure 10. — Ten-bay truss (members are circled, nodes are not). 


TABLE IX — SIMPLY SUPPORTED TRUSS 
DESIGN RESULTS 


Results 

MFUD 

FUD 

SUMT 

FD 

Optimum weight, lb 

734.15 

808.74 

719.69 

782.52 

Member area, in . 2 





a 2 

2.38 

2.73 

2.54 

2.68 

A„ 

3.46 

4.10 

3.33 

4.16 

A 25 

5.46 

6.29 

5.73 

5.23 

A 35 

5.23 

5.98 

5.03 

4.99 

A,, 

1.00 

1.23 

1.00 

1.93 

Active constraints 





Stress 

3 

— 

13 

11 

Displacement 

2 

2 

2 

2 
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Example 9: Sixty-Bar Trussed Ring 

A ring idealized by 60 bar members (ref. 23) subjected to 3 
loads is depicted in figure 1 1 . It has 60 stress and 6 displacement 
constraints for each load case. The 60 bars’ areas were linked 
to obtain 25 independent design variables. Table X presents the 
optimum designs obtained for the ring. For this example, 
MFUD, SUMT, and FD results were in good agreement, and 
the active constraints for each method numbered 19, 12, and 15, 
respectively. 


4 



10 


Figure 1 1 . — Sixty-bar trussed ring (members are circled, nodes are not; 
R q = outer radius; R f = inner radius). 


TABLE X.— SIXTY-BAR TRUSSED RING 
DESIGN RESULTS 


Results 

MFUD 

FUD 

SUMT 

FD 

Optimum weight, lb 

308.07 

324.23 

308.96 

308.93 

Member area, in. 2 





a 5 

0.59 

0.63 

038 

0.57 

Ajq 

1.84 

1.61 

1.94 

1.85 

A 15 

0.77 

0.57 

0.69 

0.71 

A 20 

0.97 

1.04 

1.07 

1.07 

A25 

1.16 

1.34 

1.15 

1.15 

Active constraints 





Stress 

18 

— 

11 

14 

Displacements 

1 

1 

1 

1 
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Example 10: Geodesic Dome 

An aluminum geodesic dome (ref. 17) idealized by 132 bar 
members and subjected to a single load is depicted in figure 1 2. 
The areas of the bars in the dome were linked to obtain the seven 
independent design variables. The designs obtained for 132 
stress and 1 displacement constraint are given in table XL For 
this example, the MFUD weight of 1 19.44 lb was lighter than 
the FD weight by 5.52 lb, but heavier than the SUMT weight by 
0.8 lb. The traditional FUD weight was 560 percent heavier. 


.un ji tun uhl 

j\ a A 

91 (92 93 94 95 96 97 S 


-^- 73 > ^ -74 “^75 ^ '! 
JO 43 (44 ; 45 : 46; 47 48 76 


/ > / N , / V / \ / X/ 30,0 in. 

129 89) ( 71 .) 42 13 14 15 16 33 50 77 1 02' V/, \ 

1 28-^- 70 30 ~~X~~ 7 (1 7 ^ (51 ) ^ 1 03“ 

127 88 69 41 29 i (12) (l ) \ 2j (8 ) 18 34 52 78 1 04 

1 26 - ^ "68 ^ -28 6 3 “^“19 ^ 53 

125 87 67 40 27 1 1 ) 15 4 ( 9 20 35 J (54 79 106 

■\66 ) 26 -^“1 ^ 55 ^ -107- 

123 86 65 39 25 24 23 (22; 36) (56 80 (108 

1 22 ^ 64 ^ -(38; ^ ( 37 r 

121 85 63) (62 61 60 59 58 (81 ) 1l0 

1 20-^- 84 83 -^82 1 1- 

119 118 117) 116 115 114 113 112 

V/iMl/rUiM TTT 


240.0 in. 



Figure 12. — Geodesic dome (members are circled, nodes are not). 


TABLE XI — GEODESIC DOME DESIGN RESULTS 


Results 

MFUD 

FUD 

SUMT 

FD 

Optimum weight, lb 

119.44 

676.01 

118.65 

124.% 

Member area, in. 2 





A, 

0.52 

1.68 

0.55 

0.90 

A, 

A, 

0.31 

1.70 

0.29 

0.29 

0.29 

1.74 

0.29 

0.29 

A, 

0.30 

1.76 

0.30 

0.30 

Active constraints 





Stress 

46 

— 

46 

50 

Displacements 

1 

1 

1 

1 
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TABLE XII —SUMMARY OF NORMALIZED WEIGHTS 


Problem 


Normalized weight 

under — 



Stress and displacement 
constraints 

Stress con 
onl; 

straints 

i 

Optimum 3 

MFUD 

FUD 

Optimum 

FSD 

Three-bar truss 

1.0 

1.00 

1.22 

1.0 

1.0 

Five-bar truss 



1.00 

1.38 

- 

- 

- 

- 

Tapered five-bar truss 



1.00 

1.00 

1.0 

1 . 

0 

Ten-bar truss 



1.02 

1.14 





Tapered ten-bar truss 



1.00 

1.03 





Twenty-five-bar truss 



1.00 

1.06 





Simply supported truss 



1.02 

1.12 





Sixty-bar trussed ring 



1.00 

1.05 





Geodesic dome 

1 


1.01 

5.70 

1 





“Normalized with respect to optimum weight obtained by SUMT. 


Step Length and Gradient 

Gradient and step length concepts are used differently in the 
MFUD method and in an optimization method. In the MFUD 
method, step length is a reduction factor and is assigned a value 
such as t = 0.5. The factor guards against rapid change in the 
updated design for which analysis may no longer be valid. Step 
length in optimization is determined from a one-dimensional 
search for a profile that is contrived or assumed by using local 
information, including the gradient vector. The step length 
reduction factor in the MFUD and the step length determination 
in optimization are quite different. In the MFUD, the gradient 
information is used to separate critical design variables. This 
separation — and consequently gradient calculation — needs to 
be carried out a few times for the entire design process: for 
example, at initial design, at final design, and at some interme- 
diate iterations. In optimization, gradient information is used to 
generate a search direction and subsequent calculations. 

Discussion 

Table XII summarizes the normalized weights obtained by 
SUMT, MFUD, and FUD for all nine examples with stress and 
displacement constraints. Solutions for stress constraints only 
were generated by using optimization and FSD methods for all 
the examples and are also included in table XII. From these 
results the following observations can be made: 

(1) For stress constraints only, the FSD’s generated by the 
stress-ratio technique are identical to the optimum designs 
obtained with mathematical programming techniques. 


(2) For both stress and displacement constraints, designs 
generated by the MFUD method are in close agreement with the 
optimum designs — only 1 or 2 percent variation (which can 
largely be attributed to the values of the convergence parameter 
of the optimization algorithms). 

(3) For stress and displacements constraints, the traditional 
FUD method typically produced overdesigns as expected. 

Concluding Remarks 

A modified fully utilized design (MFUD) method has been 
developed for the design of structures with both stress and 
displacement limitations. In the development of the MFUD 
method, the Integrated Force Method was found to be the 
appropriate analysis tool. The MFUD method has been verified 
through successful solutions of a number of design examples. 
It alleviates the overdesign limitation associated with the tradi- 
tional fully utilized design method. The MFUD method has the 
potential to become an industrial design tool for practicing 
engineers, since this simple approach can generate 
designs comparable to those produced with design optimiza- 
tion methods based on difficult nonlinear mathematical 
programming techniques. A fully utilized design, which by 
definition is a design wherein the number of active constraints 
equal or exceed the number of design variables, represents the 
optimum condition. The MFUD method needs to be developed 
for dynamic constraints and nontruss type structures. 


12 


NASA TM— 4743 




Appendix A 

Analysis Tool for the Modified Fully Utilized Design 


Either the displacement method or the force method can be 
used to develop the modified fully utilized design (MFUD) 
technique. The IFM (integrated force method) offers certain 
advantages in the development of MFUD. In this appendix, the 
basic IFM equations are summarized, and the suitability of the 
IFM for design is explored. 

Equations of the Integrated Force Method 

In the IFM, the n internal forces {F} are considered the 
primary unknowns, and these can be obtained from the solution 
of its governing equation as follows: 

[S]{F} = {P*} (20) 

where [S] is the n x n governing matrix, and {P*} is the 
n-component thermomechanical load vector. 

The m displacement components {X} can be obtained from 
the forces by back-calculation as follows: 

{*} = [J]{[Gl{F} + {j8} 0 } (21) 

where [J] is the m x n deformation coefficient matrix ([/] = m 
rows of [S]' T ), [G] is the n x n flexibility matrix, and {/3} 0 is the 
n component for the initial deformation vector. 

Notice that the IFM provides two sets of equations (eqs. (20) 
and (21)), one for the calculation of forces and another for 
displacement computations. For more details on IFM, see 
references 14 to 16. 

Equations of the Displacement Method 

The displacement vector (X) of dimension m is the primary 
unknown in the displacement method, and it can be obtained 
from its governing equation 


Suitability of Analysis Methods for Design 

The suitability of the force and the displacement methods for 
the development of MFUD can be illustrated by considering the 
example of a three-bar truss (fig. 2). The IFM governing 
equations for forces for the three-bar truss have the following 
explicit form: 



where F, is the force in the ith member; f. = {€^E f A i ) is the 
flexibi lity of the ith member for Y oung’ s modulus F-, length F-, 
and area Af 9 and P x and P are the applied load components. 

The two displacement components (X ; , for i = 1 , 2) of the 
three-bar truss can be obtained from the forces by using the 
following equation: 

yftf 1 + _1_) + V2 f 1 | 1 ) 

j X] 1 _j_ f\\h h) hyh h) ft l/i h) 
|x 2 J y __%/2_ V2 ( i | 1 ) V2 
hh h A ft) Aft 

A 0 0 1 [ F\ 

x \ ° f 2 0 U F 2 \ (24) 

0 0 / 3 J[f 3 ^ 


[*]{*} = {/>} (22) where 


where [K] is an m x m stiffness matrix, and {P} is an 
m- component load vector. 

Unlike the force method, the displacement method does not 
provide two sets of equations, one for displacements and 
another for forces. Instead, from nodal displacements a series 
of operations (such as determining the field displacement func- 
tion, computing the strain by differentiation, and then calculat- 
ing the stress by using Hooke’ s law) are carried out to determine 
internal forces. 


v = — — + 2 — — + — - — (25) 

/ 1/2 Aft ft ft 

Consider the design of the truss for stress limitations o oi , for 
i - 1 , 2, 3, and displacement limitations X oi , for / = 1 , 2. The two 
IFM equations (eqs. (24) and (25)) can be written in terms of 
member areas (A^ A 2 , A 3 ) as follows. 
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For stress limitations: 


g oi 

V2 

q 01 

42 

a oi/i 

42 


0 

~°02 

~ a (nf2 


q 02 

V2 

a 03 

42 

a Q3/3 

V2 


'V 


Px 

Aj 

> = - 

p y 

A 2 


0 


For displacement limitations: 


Even for displacement constraints only, the three areas cannot 
be directly determined from the solution of two stiffness 
equations. The stiffness formulation is not the most appropriate 
analysis tool to develop a direct design formulation. 

Displacement Sensitivities 

The n x m displacement sensitivity matrix [VX] required in 
the MFUD method can be obtained in explicit form with the 
IFM (ref. 19) as 


A 


1 El 


r 42 ' 


[VX] = [[y][C][^]-[y][G][q3]] T (29) 





V X 2 J 


where 


(27) 


m= Diag{F x ,F 2 F n 


A e 3 Fs ( -42 ' 

3 £3 {x x + x 2 j 

In equation (27), if a displacement component X i exceeds its 
allowable value X Q •, then that component should be replaced by 
the allowable value. 

An iterative solution of equations (26) and (27) can provide 
a design for the three-bar truss that accommodates both stress 
and displacement constraints. However, convergence diffi- 
culty can occur if this solution is used for more general trusses. 
The modified fully utilized design method developed in this 
paper is more suitable for such applications. 

The displacement method, on the other hand, provides only 
two equations in terms of the three bar areas: 


2 


E\A\ 

E 


EjA 3 

E-iA 


.3/13 


<1 

E X A X 


'1 


E\A\ E 3 A 3 

t ') £ 3 

+ 2 - — — ^ — — 


E2A2 E^A 


3 /1 3 




(28) 


[C] = Diag. 


[G]=Diag\ 


z4i zh z4n 

E X A^' E 2 Al""' E n A 2 n . 


*1 ^2 


(30) 


E \A\ E 2 A 2 E n A n 


m=[s\ 


,-1 


[0] 

_[C][C] 


m 


In the previous equation, Diag . designates a diagonal matrix, 
and [C] is an (m - n) x n compatibility matrix of IFM. The 
first term in equation (29) accounts for changes in member 
flexibility, whereas the second term accounts for the changes in 
member forces with respect to member areas. However, Berke 
(ref. 9) has shown that the second term is identically equal 
to zero, which has also been numerically verified. The first term 
in equation (29) is equivalent to equation (12), which is used to 
develop the MFUD method. 
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Appendix B 

Optimality of a Fully Stressed Design 


A fully stressed state is reached when all members of a truss 
are utilized to their full strength capabilities. Historically, such 
a design was considered optimum, but recently this optimality 
has been questioned because the weight of the structure is not 
used in the design calculations. This appendix examines 
optimality of the fully stressed design (FSD) with analytical 
and graphical illustrations. Solutions for a set of examples 
obtained by using FSD and optimization methods confirm the 
optimality of FSD. FSD, which can be obtained with little 
calculation, can be extended to displacement constraints and to 
nontruss-type structures. 


The optimum solution — variables (A° pt for / = 1, 2, ... , /t), 
minimum weight ( W opt ), and active constraints (gj a - 0, 
j = 1 , 2, . . . , n ) — can be obtained by using one of several 
optimization methods (see refs. 13 and 29). In optimization 
methods, both the weight function and the constraints partici- 
pate. In FSD, only the constraints are solved iteratively to 
obtain the design variables, without any reference to weight. 
The FSD weight (Wf sd ) is back-calculated from the areas. That 
FSD need not be optimum (i.e., A^ sd f. A° pt for i = 1 , 2, . . . , n, 
and Y/f sd t- W opt ) is a popular misconception. 


Introduction 

Researchers are baffled by two conspicuous attributes of 
FSD: the good numerical results obtainable with FSD; and the 
merit function, or weight function, of the structure, which is not 
taken into consideration. Optimization proponents think that 
FSD need not represent the optimum since the good FSD results 
are considered simply special cases. Practicing engineers 
believe that when all the members of a truss (or structure) are 
utilized to their full strength capabilities the design can no 
longer be improved. They, however, cannot offer a mathemati- 
cal proof supporting the optimality of FSD. This dilemma has 
persisted since the sixties (refs. 1, and 24 to 28). Here, an 
attempt is made to alleviate the confusion. The optimality of 
FSD is examined in four sections: the problem is defined; 
optimality is discussed; numerical examples follow; and dis- 
cussions and a summary are presented. 


Truss Design Problem 

Consider an n-bar truss with n member areas as design 
variables subjected to q load conditions. A fully stressed state 
(of FSD) is reached when each members’ stress equals allow- 
able strength o Q . This design can be cast as the following 
mathematical programming problem: Find n variables A i for 
i = 1 , 2, . . . , n to minimize weight w = ^ p,l,A subjected to 
nq stress constraints 1=1 


8i 


°0 


i = 1, 2 , nq 


(31) 


Optimality of the Fully Stressed Design 

The Lagrangian functional obtained by adjoining the active 
constraints to the weight function is used to examine the 
optimality of FSD: 

3!(M.M) = w (M}) + X (32) 

active set 

where (*) indicates the active constraints and {^} the multi- 
pliers. The variables and the multipliers can be obtained from 
its stationary condition: 

Vw({A}) + £ W({A}) = {0} (33) 

active set 

£i({^}) = {0} [gi within the active set j (34) 

Equations (33) and (34) yield the optimum solution. 

The optimality of FSD is considered by examining three rela- 
tions between the design variables and the active constraints. 

Case 1: There are more active constraints than design 
variables. 

Case 2: There are an equal number of active constraints and 
design variables. 

Case 3: There are fewer active constraints than design 
variables. 

The three-bar truss (fig. 2) subjected to two load conditions, 
with three design variables, six stress constraints, and weight as 
the merit function, is used for illustration. 
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Case 1: More Active Constraints Than Design Variables 

Geometrical solution . — Consider an optimum solution with 
n variables and (n + v) active constraints. The optimal solution 
is at the intersection of any n out of the ( n + v) active constraints. 
The remaining v are follower constraints passing through the 
optimal point. For the truss with three design variables, assume 
an optimal design with four active constraints, g v g 3 , g 5 , and g 6 
(fig. 13). Three constraints (g 3 , g 5 , and g 6 ) are sufficient to 
establish the optimal point. The follower constraint (gj) can be 
neglected without any consequence. From a geometrical con- 
sideration, the inclusion of a maximum of n active constraints 
is sufficient to establish the optimal design. The weight func- 
tion is not essential when v > 0. 

Analytical solution . — The (2 n+v) unknowns (being n vari- 
ables and (n + v) multipliers) can be determined as the solution 



Figure 1 3. — Three active constraints (sufficient to 
determine optional point) and two follower 
constraints. 

to equations (33) and (34). An uncoupled strategy is to solve for 
the n design variables from any n of (n + v) constraint functions 
given by equation (34). Values for other variables and the 
weight function can be back-calculated. Summarizing, when 
active constraints exceed design variables, the optimum can be 
obtained from the solution of a set of n active constraints. 

Case 2: An Equal Number of Active Constraints and 
Design Variables 

An optimal solution with n variables and n active constraints, 
by definition, represents a fully stressed design. The stationary 
condition of the Lagrangian (eqs. (33) and (34)) represent 2n 
equations in 2n unknowns. The uncoupled equation (34), being 
n constraint equations, can be solved for the n design variables. 
The n multipliers and optimum weight can be back-calculated. 
For the truss, the solution of three constraints will yield the 
design variables. The optimum weight and the multipliers can 
be back-calculated from equations (31) and (34) respectively. 


When the number of active constraints equal or exceed the 
number of design variables, the solution of the active con- 
straints (i.e., eq. (34)) provides the design variables. The design 
thus obtained is both fully stressed and optimum. 

Case 3: Fewer Active Constraints Than Design Variables 

An optimum solution with fewer active constraints than 
design variables is not a fully stressed design. For the three-bar 
truss, assume two active constraints (g l and g 2 ) given by equa- 
tion (34). The two constraint equations are expressed in terms 
of three unknown design variables. Although equation (34) is 
independent of Lagrangian multipliers, it does not have suffi- 
cient quantity for a solution of the three design variables. Thus, 
both equations (33) and (34), which are coupled in variables, 
multipliers, and weight gradient, must be solved simulta- 
neously to generate the optimum solution. The gradient of the 
weight function and the multipliers are required to calculate the 
design variables. In other words, only when the number of 
active constraints is fewer than the number of design variables 
do both the constraints and the weight function participate. 
Mathematical programming methods address this situation in 
particular. Practical truss design, however, more frequently falls 
under Cases 1 and 2. 

Design of a Truss Under a Single Load Condition 

For an indeterminate truss under a single load condition, a 
full stress state may not be achievable because of the compat- 
ibility condition (refs. 27, 28, and 30). Take, for example, an 
n-bar truss with r redundant members. If its FSD is attempted 
without restricting the lower bound of the member areas, then 
the design will degenerate to a determinate structure that, of 
course, will be fully stressed and optimum. If, however, a mini- 
mum bound A min is specified for member areas, the resulting 
design will have (n - r) fully stressed members with ( n - r ) 
active stress constraints and r member areas that reach the 
minimum bounds of A min . These properties, from an analytical 
viewpoint, become equivalent to n active constraints consisting 
of (n - r) stress constraints and r lower bound side constraints. 
Since there are n design variables, this example falls under 
Case 2. In other words, the design of a truss under a single load 
also represents the optimum design. 

A fully stressed design state can be defined in terms of two 
indices, Index stress and Index 311 : 


. , stress 

Index 


(number of active stress constraints) 
(number of independent design variables) 


^ (number of active stress constraints + number of active bounds) 
Index = 


(number of independent design variables) 
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Index - maximum |lndex stress , Index 111 j 

For analytical purposes, a fully stressed state is reached when 
the Index > 1 . 


Numerical Examples 

Examples are separated into a first example and a group of 
problems. The first example, with several subcases, examines 
the role of the weight function when the number of active con- 
straints exceed or equal the number of design variables (as in 
Cases 1 and 2). The second group of examples compares stress- 
ratio-based FSD’s with their optimum designs. Two optimizers, 
SUMT (Sequential Unconstrained Minimization Technique) 
and IMSL (i.e., the Sequence of Quadratic Programming tech- 
nique of IMSL) are used to solve examples in group 2. The 
figures and descriptions for the examples are not given here, but 
can be found in references 13 and 29. All the examples were 
solved in a controlled environment on the NASA Lewis Cray 
Y-MP computer. 


Example 1 

A three-bar truss (see fig. 2) is used to illustrate that the 
weight function does not influence the optimum design when 
the number of active constraints equal or exceed the number of 
design variables. The truss is subjected to two load conditions 
and has a total of six stress constraints, three per load condition. 
The optimum solution for an aluminum truss with equal weight 
densities of 0.1 lb/in. 3 for its three bars was obtained by using 
several optimization algorithms. The optimum solution for the 
problem is optimum weight W° pt = 133 lb; design variable 
A opt - (3.29, 3.99, 3.32) in. 2 ; and four active stress constraints 

C?l' #5> and Rf)- 

Fully stressed design . — The stress-ratio-based FSD pro- 
duced the optimum design. The weight coefficients were 
changed over a wide range, from 0.1 lb/in. 3 for aluminum to 
300 lb/in. 3 for a fictitious material. The design and active con- 
straints obtained by FSD remained the same since the weight 
does not participate in the calculations. The FSD weight, how- 
ever, was back-calculated; it is shown in table XIII. 

SUMT optimizer . — Solutions for five different weights were 
attempted by SUMT (see table XIII). The SUMT optimizer 


TABLE XIII.— OPTIMUM DESIGNS OF THREE-BAR TRUSS WITH 
DIFFERENT MATERIALS 


[Number of active constraints exceed number of design variables.] 


Method 

Cost 

coefficients 

Member areas 

Active 

constraints 

Optimum 
weight, lb 

FSD 

0.1 

0.1 

0.1 

3.30 1 

3.99 1 

3.32 1 

81.83 ’85 *86 

1.33x10* 


3 

6 

8 









7.53x10* 


6 

12 

18 









1.60x10* 


16 

13 

25 









2.43x10* 


1 

200 

300 









2.20x10* 

SUMT 

0.1 

0.1 

0.1 

3.291 

3.986 

3.323 

81,83*85*86 I 

1.33 xlO 2 


3 

6 

8 

3.299 

3.998 

3.299 



7.53x10* 


6 

12 

18 

3.299 

3.997 

3.298 



1.59x10* 


16 

13 

25 

3.298 

3.998 

3.299 



2.43x10* 


1 

200 

300 

67.068 

9.111 

0.001 



1.92x10 s 

ISML 

0.1 

0.1 

0.1 

1.000 

1.000 

1.000 

(a) 

3.83x10* 


3 

6 

8 

3.299 

4.000 

3.299 

8| *83 ’85 *86 

7.53x10* 


6 

12 

18 









1.60x10* 


16 

13 

25 









2.43x10* 


1 

200 

300 









230x10* 

OPTM1 

0.1 

0.1 

0.1 

3.313 

3.971 

3.323 

81*83 *85*86 

1.33x10* 


3 

6 

8 

3.309 

3.963 

3.334 



7.55x10* 


6 

12 

18 

3.309 

3.962 

3.335 



1.60x10* 


16 

13 

25 

3.308 

3.961 

3.336 



244x10* 


1 

200 

300 

3.300 

3.967 

3.328 



221x10 s 

SQP 

0.1 

0.1 

0.1 

2.335 

2.503 

2.505 

(a) 

9.35x10* 


3 

6 

8 

2.334 







5.32x10* 


6 

12 

18 

2.334 







1.14x10* 


16 

13 

25 

2.335 







1.74x10* 


1 

200 

300 

2.335 







1.57 x10 s 

OPTM2 

0.1 

0.1 

0.1 

3.199 

2.556 

5.102 

81*85 

1.42x10* 


3 

6 

8 

3.501 

2.402 

4.682 

(a) 

8.22x10* 


6 

12 

18 

3.635 

2.353 

4.563 

86 

1.75x10* 


16 

13 

25 

3.65 7 

2.316 

4.467 

86 

271x10* 


1 

200 

300 

3.684 

2.336 

4.512 

86 

2.39x10 s 


a No active constraints. 
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TABLE XIV — PERFORMANCE OF FULLY STRESSED DESIGN VERSUS OPTIMIZATION METHODS 


Problem 

Load 

condition 

Independent 

design 

variables 

Design 

method 

Normalized 

weight 

Number of 
active stress 
constraints 

Number of 
active side 
constraints 

Index 111 

Index strcss 

Normalized 
CPU time 

BAR3 

1 

3 

FSD 

1.000 

2 

1 

1.0 

0.667 

1.000 




IMSL 

1.000 

2 

1 

1.0 

.667 

1.596 




SUMT 

1.000 

2 

1 

1.0 

.667 

6.263 

BAR5 

2 

5 

FSD 

1.000 

7 

0 

1.4 

1.400 

1.000 




IMSL 

1.000 

7 

0 

1.4 

1.400 

.527 




SUMT 

1.000 

8 

0 

1.6 

1.600 

1.288 

BAR10 

1 

10 

FSD 

1.000 

6 

4 

1.0 

0.600 

1.000 




IMSL 

1.000 

6 

4 

1.0 

.600 

1.362 




SUMT 

1.001 

6 

4 

1.0 

.600 

2.859 

BAR25 

2 

8 

FSD 

1.000 

11 

3 

1.75 

1.375 

1.000 




IMSL 

1.000 

11 

3 

1.75 

1.375 

8.688 




SUMT 

1.001 

11 

3 

1.75 

1.375 

14.213 

DOME 

1 

12 

FSD 

1.000 

188 

5 

16.08 

15.667 

1.000 




IMSL 

.983 

192 

5 

16.42 

16.000 

.743 




SUMT 

.984 

192 

4 

16.33 

16.000 

.740 

RING_A 

3 

25 

FSD 

1.000 

40 

0 

1.60 

1.600 

1.000 




IMSL 

.999 

38 

0 

1.52 

1.520 

5.476 




SUMT 

1.000 

38 

0 

1.52 

1.520 

13.101 

RING_B 

1 

60 

FSD 

1.000 

52 

16 

1.133 

0.867 

1.000 




IMSL 

1.000 

52 

16 

1.133 

.867 

2.882 




SUMT 

1.003 

52 

8 

1.0 

.867 

5.569 

TOWER.A 

1 

252 

FSD 

1.000 

117 

135 

1.0 

0.464 

1.000 




EMSL 

.999 

117 

131 

.984 

.520 

57.249 




SUMT 

1.000 

117 

139 

1.016 

.551 

81.442 

TOWER_B 

2 

252 

FSD 

1.000 

165 

97 

1.040 

0.655 

1.000 




IMSL 

1.000 

165 

98 

1.044 

.655 

48.031 




SUMT 

1.000 

165 

99 

1.048 

.655 

t£ 59.557 


converged to the optimum solution for the first four cases, pro- 
ducing the correct optimum weight and an identical set of four 
active constraints. For the fifth case, SUMT converged to an 
eccentric local optimum design with two, instead of four, active 
constraints. 

IMSL optimizer . — This optimizer produced correct solu- 
tions for the last four cases. For the first case, no active con- 
straints were produced and the IMSL solution was unsatisfactory. 

Table XIII also shows solutions obtained with other optimi- 
zation methods. The FSD method provided successful solu- 
tions for all five cases. The success rate for optimization 
methods varied. 

Example Set 2 

Nine examples were solved by using FSD along with SUMT 
and IMSL optimizers. The normalized results with respect to 
the FSD answers are depicted in table XIV. 

Consider the 25-bar truss, referred to as B AR25 in table XIV. 
It is subjected to two load conditions. Its 25 areas are linked to 
obtain eight independent variables. Since, at optimum, 11 
stress constraints are active, a fully stressed state has been 
reached. The FSD, IMSL, and SUMT methods produced iden- 
tical optimum solutions for the example with different CPU 
times. Optimizers IMSL and SUMT were, respectively, 1 0 and 
18.5 times more expensive than the FSD method. Solutions to 
the other eight problems followed the pattern of the 25-bar 
truss, with minor variations. 


Discussion 

For a truss, if a fully stressed state can be reached (i.e., the 
number of active constraints exceed the number of design 
variables), then such a design can be handled satisfactorily with 
the stress-ratio-based FSD method. Optimization techniques 
for such problems can be computationally expensive and un- 
necessary. 

In special circumstances a practical structural design may be 
associated with fewer active constraints than design variables. 
Such a design is likely to represent an overdesign condition, 
which can be alleviated by relaxing some of the nonactive 
constraints. If, however, there are fewer active constraints than 
design variables, then the design is not fully stressed; here, non- 
linear programming optimization methods can be useful. For 
such problems the stress-ratio-based design can differ from the 
optimum design, especially when weight densities for truss 
members are different. 

When the fully stressed design is extended to include dis- 
placement constraints, it is called a fully utilized design (FUD). 
FUD, which can produce overdesign conditions, has been 
modified to give a method that produces a satisfactory design 
for stress and displacement constraints (ref. 31). 

The FUD method has been extended in reference 32 to 
nontruss-type structures. 
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Summary 

A fully stressed design is optimum when a full stress state can 
be achieved. At optimum, when the number of active con- 
straints equal or exceed the number of design variables, then 
such a design can be obtained by simply using a stress-ratio 
algorithm without any consideration to the weight function. 
The stress ratio algorithm can produce a fully stressed design in 
a small fraction of the calculation time required by the design 
optimization methods. The fully stressed design method may 
have the potential for extension to nontruss-type structures and 
nonstress constraints. 
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